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ABSTRACT 

We analyze the relationship between maximum cluster mass, and surface densities of total gas 
(Egas), molecular gas (Shs), neutral gas (Shi) and star formation rate (Ssfr) in the grand design 
galaxy M 51, using published gas data and a catalog of masses, ages, and reddenings of more 
than 1800 star clusters in its disk, of which 223 are above the cluster mass distribution function 
completeness limit. By comparing the 2-D distribution of cluster masses and gas surface densities, 
we find for clusters older than 25 Myr that Afsrd oc E^^^*^'^, where Afsrd is the median of the 
5 most massive clusters. There is no correlation with Sgas, Sh2, or Ssfr- For clusters younger 
than 10 Myr, Ahrd oc E^f^^'^, M^rd oc Sgag^*^'^; there is no correlation with either or Ssfr- 
The results could hardly be more different than those found for clusters younger than 25 Myr 
in M 33. For the fiocculent galaxy M 33, there is no correlation between maximum cluster mass 
and neutral gas, but we have determined Ma^d oc S^^s'^''-^; Msi-d oc Sna"^"'^; oc Ssfr°'^- 

For the older sample in M 51, the lack of tight correlations is probably due to the combination 
of the strong azimuthal variations in the surface densities of gas and star formation rate, and the 
cluster ages. These two facts mean that neither the azimuthal average of the surface densities at 
a given radius, nor the surface densities at the present-day location of a stellar cluster represent 
the true surface densities at the place and time of cluster formation. In the case of the younger 
sample, even if the clusters have not yet traveled too far from their birthsites, the poor resolution 
of the radio data compared to the physical sizes of the clusters results in measured Es that are 
likely quite diluted compared to the actual densities relevant for the formation of the clusters. 



Subject headings: galaxies: star clusters 
galaxies: individual (M 51, NGC 5194) 



galaxies: ISM — galaxies: spirals — stars: formation 



1. Introduction 

This is the second in a series of papers devoted 
to investigate whether the maximum young clus- 
ter masses observed in a galaxy reflect physical 
processes that determine them or if, instead, they 
constitute an effect of random sampling statistics. 

A correlation between gas surface density and 
embedded maximum cluster mass (i.e., Med, max oc 
Egas) has been proposed as an ansatz to explain 
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the existence of the Ha cut-off in disk galax- 
ies, while at the same time, as can be inferred 
from UV da ta, there is no actual star forma- 
tion cut-off (IPflamm-Altenburg fc Kroupal liooj 
Pflamm- Altenburg et al.l l2009l) . Similarlv. a cor- 



rel ation Med. ma x oc Egpj^ has been pr oposed 
by iBillett et all (|2002l) and iLarsenl (|2002l ). aris- 
ing from the st ar formation law (Esfr oc Eg^^ 
IKennicutt 1998() . in conjunction with the assump- 
tion of pressure equilibrium between the ambient 
ISM and the cluster-forming cloud cores. In this 
case, 2/3 < rj < 2, with the lower value expected 
for equal (volume) density cluste rs, and the h igher 
one for clusters with equal size ( LarsenI 2002 ) . 



In the first paper of the series (jGonzalez-Lopezlira et al 
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2012I ). we analyzed radially averaged distributions 
of gas surface densities and of cluster masses in 
the flocculent galaxy M 33. Here, we study the 
grand-design galaxy M 51. Star cluster data are 
presented in Section 1, and ISM data in Section 2; 
both cluster and gas data are taken from the lit- 
erature. The results and a comparison with M 33 
are discussed in Section 3. Section 4 contains our 
summary and conclusions. 

2. Star cluster data 

We rely on two catalogs of star clusters brighter 
than Vp555Vi/ = 23 mag, detected by Hwan g &: Le'3 
(|2008l ). based on Hubble Space Telescope (HST) 
Advanced Camera for Surveys (ACS) observa- 
tions acquired by the Hubble Heritage Team. 
Hwang fc Led compiled two catalogs, one of 2224 
"Class 1" clusters, that appear circular and well 
isolated, and one with 1388 "Class 2" clusters, 
which have elongated or irregular shapes and/or 
close neighbors. These original catalogs list, 
among other properties, cluster RA and DEC, 
half-light radius R^s, (B - V) and (V - I) col- 
ors, aperture corrected V mag, and photometric 
errors, all in the HST equivalent of the Johnson 
system. From the r ecovery of artificial clusters, 
Hwang fc Led ( 2008 ) determined that the com- 
pleteness of clusters with V < 22, mag is higher 
than 80%, regardless of the clustering properties 
of the field. 

Subsequently, iHwang fc Le^ (|2010l) calculated 
masses, ages, and reddening for the clusters, via 
the comparison of their broad ban4ll spectral en- 
ergy distributio ns (SEDs) with theoretica l SEDs 
derived from the Bruzual fc Chariot (l2003h stellar 
population synthesis models. The models used by 
Hwang fc Led ( 2010l ) have a Salpeter initial mass 



function (IMF), with lower and upper mass cut- 
offs of 0.1 Mq and 100 Mq, respectively; ages be- 
tween 1 Myr and 15 Gyr; and solar metallicity. 
Each one of the 182 model spectra was reddened 
by E{B — V) values between 0.0 and 0.6 in 40 
0.015 mag stepsH The age and reddening of each 
cluster were determined simultaneously, from a 



1 The Hubble Heritage ACS data set for M 51 com- 
prises images in the FAZhW , FbbbW , and FSliW filters. 
iHwang &: Led l(2010.} complemented these with archival 
F336V1/ Wide Field Planetary Camera 2 (WFPC2) data. 

■^This range in color excess includes the values obtained for 
M 51 clusters by previous studies; most clusters have E(B— 



fit of the reddened model SEDs to the photomet- 
ric data. The original catalogs were thus reduced 
to 1125 Class 1 and 835 Class 2 objects with ac- 
ceptable fits. For each cluster, a mass was finally 
assigned from the mass to light ratio MjLy of 
the best-fit model and the y-band cluster mag- 
nitude, taking the measured extinction into ac- 
count aiid_ado£tingadistance to M 51 of 8.4 
Mdc dFeldmeier et al.l[l997[ )FI The iHwang fc Lee 
( 2010l ) catalogs list, for each cluster, log age and 
errors, log mass, color excess and errors, and 
value of the best- fit. These catalogs do not provide 
an error in the mass determination. One possible 
estimate is given by the MjLy ratios at the mini- 
mum and maximum ages (the age minus and plus 
its error, respectively) of each cluster. Except for 
24 (out of 1960) clusters (half of them less mas- 
sive than 5 x IO'^Mq), this error is smaller than 
0.15 dex (35% of the mass) in all cases. The dom- 
inant source of error is probably the uncertainty 
in the distance modulus (0.6 mag, or a factor of < 
2 in mass). It would affect all the masses equally, 
though, so it would not change the exponent of 
any relationship between cluster mass and gas or 
SFR surface density. 

The age distribution of the clusters shows three 
peaks, that likely signpost bursts of enhanced star 
formation rate, respectively, at stellar ages of « 
5, 100, and 200 Myr. The two older starbursts 
are better traced by the detected Class 1 objects, 
which have masses ranging from «3xl0'^to> Ix 
IO^Mq. Conversely, detected clusters in the recent 
burst belong mostly to Class 2, with significantly 
lower masses, between w 600 and 4 x IO^Mq. 

3. ISM data 

Since M 51 is a grand-design spiral with large 
azimuthal variations of gas surface density, a two- 
dimensional approach is required, as opposed to a 
comparison of individual cluster masses with ra- 
dial gas profiles, that was adequate for M 33. 

We use the HI map of M 51 data from 
Rots et all ( 1990l ): it was originally acquired with 



the Very Large Array (VLA) and has a spatial 
resolution of 13". We derive the molecular gas 



V) < 0.1 l lLee et al.ll2005l) . 

^ At this distance, l" = 40.7 pc; R25, the galactocentric 
distance of the isophote with surface brightness in the B- 
band fiB = 25 mag D"~^, is 13.7 kpc or 5'.6. 
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surface density distr ibution from the ^^C O 2-1 
mosaic obtained by ISchuster et al. ( 2007 ) with 
the IRAM-30 m telescope; it has a resolution of 
11". Finally, to calculate the SFR surface density 
distribution we rely on the total p ower 20 cm ra- 
dio co ntinuum image published bv lFletcher et al.l 
(|201ll ): it has a resolution of 15", and was pro- 
duced from VLA C- and D-array data. 

In order to get the H2, HI, and SFR surface den- 
sities from the integrated i ntensity maps, we fo l- 
low the same procedures as lSchuster et al.l (|2007| ). 
At a given position in M 51, the neutral gas sur- 
face density Ehi = mH-/V(HI) cos z, where toh 
is the hydrogen atomic mass and i — 20° is 
the inclination angle of M 51; assuming optically 
thin emission, iV(HI) = 1.82 x 10^^ J T^^dv / {K 
km s~^)cm~^. Likewise, we obtain the molecu- 
la r gas surface density as = 2 7niiN{}l2) cos 



Schuster et al.l pOOTI ) assume a CO 2-1/1-0 in- 



tensity of 0.8, and a constant CO to H2 conver- 
sion factor, X, 1/4 that of the Milky Way (MW), 
with Xmw = 2.3 X 10^° cm-^ (K km s-^)-\ 
such that the molecular hydrogen column density 
is A^(H2) = 0.25 Xmw(1/0.8) / T„,b(C0(2-l))dw, 
where Tmb is the main beam brightness temper- 
ature. The total gas surface density is Egas = 
1.36(Sh2 + Sht), so as t o inclu de helium. Again, 



following ISchuster et al.l pOOTl ). we derive Ssfr 
from the 20 cm radio continuum data|f| The 20 
cm non-thermal radio continuum is due to syn- 
chrotron radiation by supernovae-accelerated cos- 
mic rays; it hence traces recent star formation, 
and has a strong correlation with the far-infrared 
(FIR; e.g.- lHelou et al.lll985[ )FI At the position of 
each star cluster, gas surface densities are obtained 
from the HI, H2, and 20 cm mosaics averaging, re- 
spectively, the emission over 20x200" (r^ 820x820 
pc^), 18 X 180" (« 730 X 730 pc^), and 18 x 18D"E 



*EsPR = 1.53 X lO^ g— r^V; Mepc-^Gyr-i. 

^ Ssppt can be readily derived from FIR emission, given that 
a significant amount of it is due to reprocessing by dust of 
light from young stars; FIR emission also has the advan- 
tage that it does not need to be corrected for extinction. In 
fact. 'Ro sa-Gonzalez et al] (2002) find that SFRs obtained 
at optical and UV wavelengths are often overcorrected for 
extinction, since stellar Balmer absorption contaminates 
emission lines; optical/UV and FIR/mm SFR estimators 
agree well once this effect is taken into account. 

^ The limit is imposed by the pixel size of the CO image, i.e., 
6" on the side; the smallest area over which we can average 
the flux is a subarray of 3 X 3 pixels (sides of subarrays must 



Ninety-five of the clusters in the sample have coor- 
dinates that fall outside of the borders of the radio 
images, so that only 1865 objects can be compared 
with the gas and SFR densities at their location 
in the M 51 disk. 

4. Results 



Hwang fc Le3 ()2008l ) find that the completeness 



distribution is quite similar for both Class 1 and 
Class 2 clusters. Since cluster brightness in the op- 
tical diminishes with age, the completeness limit of 
the cluster mass function depends on cluster age. 
Based on thei r photometric completeness limits, 
Hwang fc Led determine that conservative limits 
over which their mass functions are not affected 
by incompleteness are w 5 x 10'^ M© for clusters 
younger than 10^, and at least IxlO^M© for clus- 
ters between 100 and 250 Myr. 

4.1. The bursts 100 and 250 Myr ago 

We will first work with the 167 clusters in the 
catalogs with at least IO^Mq, ages between 40 
and 400 MyrQ and detected CO, HI, and 20 cm 



emissi on at their present loc a tions. IChandar et al 



(|2010t) and IChandar etllj (|201l[) find that the 
cluster mass function is approximately indepen- 
dent of age (and the age function independent 
of mass) for clust ers younger than 4 x 10® yr 



in M 83 and M 51 (|Chandar et al.ll201lD . respec 



tively. S imilar results are found for th e Magellanic 
Clouds (IParmentier fc de Griif^'2008') and the An- 
tena galaxies (Fall et al. 2009). The implication is 
that th e mass-dependent dis ruption timescale is > 
2 Gyr (IChandar et al.ll2010l) . 

We take care here to compare cluster mass 
with gas and SFR surface density always in bins 
with the same number of clusters, in order to test 
whether any change in maximum cluster mass is a 
statistical, size-of-sample effect We search for 



have an odd number of pixels), centered on the position of 
each cluster. We match the elements over which we average 
the HI and 20 cm fluxes to this 18" X 18" footprint as closely 
as is allowed by the pixel sizes of the images, respectively, 
4" X 4" (arrays with 5 pixels on the side) and 2" x 2" 
(arrays of 9 X 9 square pixels). 

^The young tale of the age distribution of the clusters in 
the 10* yr old burst reaches 25 Myr, but no cluster in this 
group that is younger than 40 Myr makes the 10^ Mq cut. 

* If clusters are drawn randomly from the same mass distri- 
bution function that declines with mass, the probability of 
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correlations between cluster mass, and gas and 
SFR surface densities of the form: 

logio ^'h = ^^logio + a^, (1) 

where Mi is the mass of the i-th most massive star 
cluster, and x stands for HI, H2, total gas or SFR. 

Figures [l] [21 |3l and |4] show log of star cluster 
mass versus, respectively, log of total gas, molec- 
ular gas, star formation rate, and neutral gas sur- 
face densities. The comparison is performed in 
bins with equal number of clusters. Bins increase 
from 2 (upper left panel) to 6 (lower middle panel) , 
and the number of clusters in each bin is indicated. 
Instead of using the maximum cluster mass in each 
bin for the comparison, the median of the five most 
massive clusters (or the third most massive clus- 
ter, Msrd) in each bin is displayed as a filled circle; 
the error bars represent the interquartile range, for 
both mass and gas/SFR surface density. Since the 
uncertainty in the measurement of any individual 
mass is typically a factor of 2-3, the median of the 
5 most massive clusters in each bin should in gen- 
eral be a more statistically robust indicator of any 
existing trend (in every bin there are always more 
than 5 clusters more massive than the complete- 
ness limit). There is no correlation of Afsrd with 
the surface densities of total gas, Sgas, molecular 
gas, Eh2, or star formation rate, Ssfr- 

Conversely, there is a hint of a correlation with 
Shi. The slopes (/3hi) and intercepts (ajjj) of the 
fits in the bins with 3 or more points, as well as 
their uncertainties, are listed in Table 1. The fits' 
average can be expressed as logio M^^d — (0.4 ± 
0.2) logio Shi + (5.2 ± 0.3). 

In order to further explore whether the intrinsic 
mass distribution of clusters may be changing as a 
function of neutral gas surface density, we perform 
a Kolmogorov-Smirnov (K-S) test on the 6 differ- 
ent subsamples presented in the bottom middle 
panel of Figure ID Figure [5] shows the cumulative 
probability distributions and Ehi of each subsam- 
ple, and the Dj^k and Pj^k values for every bin pair 
j, k are given in Table 2 (higher bin number indi- 
cates larger surface density). Assuming that sam- 
ple pairs with Pj k < 0.05 are taken from different 
distribution functions with high significance, bin 1 
{red solid line), with the smallest surface density 



picking a massive cluster increases with the size of sample. 



Shi = 6.8^2^9-^0 P^"^, seems to be different from 
bin 6 {orange long dashed-dotted line), with the 
largest surface density Shi = 20.9~!il'^M(7) pc~^ 
{Pis — 0.03). The mass distributions in bins 3, 
4, and 5 appear to lie somewhere in between. No 
hard conclusion can be drawn, however, since bin 
2 {green dotted line), with the second lowest HI 
surface density, is in fact indistinguishable from 
bin 6 (P2,6 > 0.98 for the pair). 

4.2. The recent burst 5 Myr ago 

We now analyze the 56 clusters younger than 
10^ yr with mass at least 5xlO'^Af0 and higher, 
and detected CO, HI, and 20 cm emission at their 
positions. Figures IH [3 El and |9] show log of star 
cluster mass versus, respectively, log of Sgas, Sh2, 
S^SFR, and Shi. Bins with indicated, equal, num- 
ber of clusters again increase from 2 to 6; points 
are the median of the 5 most massive clusters in 
each bin. 

As with clusters in the two older bursts, there 
is no correlation between Msrd, and Sh2 or SgpR. 
The correlation with Shi is tighter and has a 
slightly steeper slope, though, with average fit for 
all the bins with 3 or more points logio M^^^ = 
(0.6±0.1) logio Shi-F(3.5±0.1). Moreover, there 
is also a hint of a correlation of Ms^d with the 
total mass surface density Sgas, with average fit 
logio Msrd = (0.5 ± 0.2) logio 5]gas + (3.4 ± 0.2). 
Fully aware that the statistical significance might 
be questionable, we show nontheless logio of max- 
imum cluster mass in the bin, Mmax, versus Shi 
for clusters younger than 10'' yr in Figure [101 
The error in log mass for each cluster was deter- 
mined as described previously in Section [51 the er- 
ror bar in Shi shows the uncertainty in the mea- 
surement at the cluster's position. The average 
of the fits in the bins with more than 3 points 
is as tight as for M^j-d, and consistent with it: 
logio M^ax = (0.7 ± 0.1) logio S^gas + (3.6 ± 0.1). 
The parameters of the fits for the younger burst 
can be found in Table 3. 

We perform once more a K-S test on the 6 dif- 
ferent subsamples presented in the bottom mid- 
dle panel of Figure [HI The test is illustrated with 
Figure 1111 and its numerical results are given in 
Table 4. Bin 1 {red solid line), with the small- 
est surface density Shi ~ 6.612.7-^0 PC~^j seems 
to be different, with high statistical significance, 
from bin 6 {orange long dashed-dotted line), with 
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Fig. 1. — M 51, 100 and 250 Myr old bursts, logio cluster mass versus logio total gas surface density, in 
bins with equal number of clusters. Filled circles: M^^d, median of 5 most massive clusters (i.e., effectively, 

the third most massive cluster) in the bins. Number of bins incrcascis from 2 (upper left panel) to 6 (lower 
middle panel); number of clusters in each bin are indicated in the panels. The bar at the bottom of each 
panel shows surface density ranges of the bins. 



Table 1 

M 51, 100 AND 250 \lYR OLD BURSTS. FiTS TO LOOm (.1/:i,,i /-'I/.:- ) VS. LOCm Hm/iM,, pc-'^) 





iVcl 


/3hi 


am 


3 


55-56 


(0.2 ±0.1) 


(5.4 ±0.1) 


4 


41-42 


(0.2 ±0.1) 


(5.4 ±0.1) 


5 


33-34 


(0.5 ±0.2) 


(5.1 ±0.3) 


6 


27-28 


(0.7 ±0.3) 


(4.8 ± 0.3) 



Note. — Col. (1): number of bins. Col. (2): number of clusters in each bin. Col. (3): best-fit slope 
(eq. 1). Col. (4): best-fit intercept (eq. 1). 
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Fig. 2. — M 51, 100 and 250 Myr old bursts, logio cluster mass versus logio molecular gas surface density. 
Symbols as in figure [1] 



Table 2 

K-S TEST D AND P VALUES, 100 AND 250 MYR OLD BURSTS 





D P 


D P 


D P 


D P 


D P 


Bin 


2 


3 


4 


5 


6 


1 


1.53 0.0187 


0.69 0.7198 


0.97 0.3003 


0.83 0.4903 


1.44 0.0318 


2 




1.11 0.1687 


0.83 0.4903 


0.83 0.4903 


0.46 0.9823 


3 






0.83 0.4903 


0.69 0.7198 


1.04 0.2337 


4 








0.28 1.0000 


0.76 0.6100 


5 










0.73 0.6530 







Note. — D and P values for bin pairs. The cell in the intersection of a row j and a column k contains the Dj^k 
and Pj^k parameters, respectively, of the comparison between the two bins indicated in the corresponding row and 
column. If Pj^k < 0.05, the null hypothesis that the clusters in the two bins are taken from the same mass distribution 
function is rejected. 
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Fig. 3. — M 51, 100 and 250 Myr old bursts, logio cluster mass versus logio star formation rate surface 
density. Symbols as in figure [T] 
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Fig. 4. — M 51, 100 and 250 Myr old bursts, logio cluster mass versus logio neutral gas surface density. 
Dashed line: linear fit. Other symbols as in figure [T] 
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K-S test, 100 and 250 Myr old bursts. Cumulative probability distributions of mass for clusters 
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Fig. 9. — M 51, burst 5 Myr ago, logio cluster mass versus logio neutral gas surface density. Dashed line: 
linear fit. Other symbols as in figure [1] 



13 



5 - N=ll-12 



O 



3-1- 



O 



3 - 



-I— I— I— I— I— r- 



-- N=9-10 



I — ^ 1 



H 1 



I I I I I I I I I I I I I I 



5 - N=18-19 



I, 



4 



-H I I I I -- I h 



-I I I I I I I I I I I I I i_ 



T 



T 



-- N=28 



I h 



I I I I I I I I I I I I I I 



-- N=14 



_l I I I I I I I I I I I I L. 



T 



-I — I — I — I — I — r- 



I — i» -O i I 



I 1 — I — I— I 



_I I I I I L_ 



0.5 1 1.5 
log Sh,/MqPC-2 



0.5 1 1.5 0.5 1 1.5 
log Sh,/MqPC-2 log Sh/MqPC-z 

Fig. 10. — M 51, burst 5 Myr ago, logio cluster mass versus logio neutral gas surface density. Empty 
diamonds: maximum cluster mass; dashed line: linear fit. Other symbols as in figure [1] 
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Table 3 

M 51, BURST 5 Myr ago. Fits to logio (M/Mq) vs. logio and a'^ 





Pm "hi 


logKj (Mir, 


/M ./) vs yui/M, pc " 


3 18-19 

4 14 

5 11-12 

6 9-10 


(0.7 ±0.1) (3.5 ±0.1) 
(0.7 ±0.1) (3.5 ±0.1) 
(0.6 ±0.5) (3.5 ±0.6) 
(0.7 ±0.4) (3.3 ±0.4) 


logio (Msrd 


/Me) vs Egas /Mq pc"^ 


3 18-19 

4 14 

5 11-12 

6 9-10 


(0.6 ±0.5) (3.3 ±0.8) 
(0.6 ±0.5) (3.3 ±0.8) 
(0.5 ±0.3) (3.4 ±0.4) 
(0.2 ±0.3) (3.7 ±0.4) 


logio (Mmax/Afo) vs Ehi/Mq pc-^ 


3 18-19 

4 14 

5 11-12 

6 9-10 


(0.8 ±0.1) (3.6 ±0.1) 
(0.8 ±0.3) (3.5 ±0.3) 
(0.7 ± 0.3) (3.6 ± 0.3) 
(1.0 ± 0.4) (3.3 ± 0.4) 



Note. — Col. (1): number of bins. Col. (2): number of clusters in each bin. Col. (3): best-fit slope 
(eq. 1). Col. (4): best-fit intercept (eq. 1). 



Table 4 

K-S test D and P values, burst 5 Myr ago 





D P 


D P 


D P 


D P 


D P 


Bin 


2 


3 


4 


5 


6 


1 


1.20 0.1108 


0.91 0.3754 


0.60 0.8646 


0.83 0.4892 


1.38 0.0438 


2 




1.12 0.1617 


1.12 0.1617 


0.81 0.5305 


0.78 0.5732 


3 






0.51 0.9575 


0.76 0.6030 


1.27 0.0778 


4 








0.51 0.9575 


1.27 0.0778 


5 










0.76 0.6030 







Note. — D and P values for bin pairs. The cell in the intersection of a row j and a column k contains the Dj^k 
and Pj^k parameters, respectively, of the comparison between the two bins indicated in the corresponding row and 
column. If Pj_fe < 0.05, the null hypothesis that the clusters in the two bins are taken from the same mass distribution 
function is rejected. 
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the largest surface density Shi — IQ.G^q ^M© 
pc~^ (-Pi, 6 = 0.04). The mass distibutions in all 
the other bins lie in between, although we notice 
that the distributions in bins 4 (cyan long dashed 
line) and 5 {magenta dashed-dotted line) are ac- 
tually closer to bin 1 than to bin 6 (Pi,4 = 0.86; 
Pi^e = 0.08; Pi^s = 0.49; Ps.e = 0.08). 

4.3. Comparison with M 33 



In a companion paper (jGonzalez-Lopezlira et al 



20121 ). we have performed a similar analysis on 
the flocculent galaxy M 33, except that we have 
compared the radial profiles of cluster mass dis- 
tribution and gas surface densities. Since M 33 
is a flocculent galaxy without a strong spiral per- 
turbation, in general one can assume that the E 
azimuthal averages provide a good estimate of 
the actual surface densities at any given point lo- 
cated at the same galactocentric distance. From 
the radial analysis, and both the average and 
median of the 5 most massive clusters in data 
bins, we derive for M 33: logio Mmcan.Srd oc 

(4.7 ± 0.4) logio ^gas; loglO Mmea„,3rd OC (1.3 ± 
0.1) logio ^Ha; loglO Mnioan,3rd OC (l.OlbO.l) log^g SsFR- 

One would also derive logio Afmoan.Srd oc (9.6 ± 
1.2) logio ^Hi, but there is a very small range 
of Shi values in the disk of M 33, as the ra- 
dial distribution of neutral gas has d T,iii/dR25 ~ 
O.2M0pc~^i?^g^ (exponential scale length = 
4.3 i?25), i.e., it is more than twice as shallow as 
in M 51, which has /ir = 2 i?25- 

Here, we present the M 33 data in a way that 
is directly comparable to our analysis of M 51. To 
this end, and since all ISM phases have roughly 
exponential profiles, we use radius as a proxy for 
surface density, and produce plots for M 33 anal- 
ogous to Figures [H El |3l and IH We note that 
the logio ^ ranges indicated by the bars at the 
bottom of the panels sometimes overlap slightly, 
owing to the departures of the surface density 
profiles from a perfect exponential; and (2) we 
omit in all cases the last bin, i.e., the one at the 
largest galactocentric dis t ance. As explained in 
Gonzalez-Lopezlira et al.l (|2012l ). this bin is dom- 
inated by the most massive cluster in the galaxy, 
whose formation conditions are badly misrepre- 
sented by the azimuthally averaged gas surface 
density. Moreover, the last bin always includes a 
radial range for which there is no gas data; these 
stop at « 0.85 i?25. The results of this exercise are 



shown in Figures [T2l[T3l[T4l and [151 and in Table 3. 
Firstly, this fresh look at the data highlights that 
a correlation does not seem to exist between M^rd 
and Shi El As for the other ISM components, we 
get: logio Mard = (3.8±0.3) logio Sgas+(0.7±0.3); 
logio Mard = (1.2 ± 0.1) logio ^h, + (3.8 ± 0.1); 
logio Mard = (0.9 ± 0.1) logio ^^sfr + (4.2 ± 0.1). 

The 2-D correlation of logio -^^3rd with total gas 
surface density, logio '^gas, has a slightly shallower 
slope than the comparison between radial profiles. 
The radial and 2-D results for molecular gas and 
star formation rate densities, on the other hand, 
are consistent within 1 a. 

5. Summary and conclusions 

We have analyzed the relationship between 
maximum cluster mass and gas surface density 
in M 51, in order to explore the suggestion that 
maximum cluster mass is determined by physical 
processes, e.g., the equilibrium pressure between 
cluster form ing cores and the ambient inte rstel- 
lar medium |Larsenll200llBillett et al.ll2002l) . and 
since the existence of such a relationship can rec- 
oncile SFR measurements derived, respectively, 
from Ha and FUV emission in galaxy disks. Fur- 
thermore, we have already found evidence in M 33 
that star formation is not a stochastic process, and 
that the range of star cluster masses is driven by 
environmental physics. 

In our work, we have us e d published ga,s data 
of M 51 (iRots et al.l [l990l: ISchuster et all 12007 : 
Fletcher et al.l 2011 ), and catalogs originally com- 
prising more than 1800 youn g star clusters in its 
disk also from the literature ( Hwang fc Leel[2008l 
2010l ). Among other information, these catalogs 
provide positions, ages, reddenings, and masses for 
the clusters. 

We have compared the 2-D distribution of the 
young star clusters with those of gas surface den- 
sity and star formation rate. To test whether any 
trend is caused by random sampling from the clus- 
ter mass function as a size of sample effect, we have 
measured the distribution of maximum mass with 
gas surface density in 2 to 6 bins with an equal 
number of clusters in each bin. 



^Notice also that the bars at the bottom of Figurc lTSl overlap 
more than a httle, a sign that Shi is far from changing 
monotonically with radius. 
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Fig. 12. — M 33, logio mass of third most massive cluster vs. logio total gas smface density. Empty triangles: 
Msrd, median of 5 most massive clusters in bin. Dashed line: linear fit. Other symbols as in Figure [T] 
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Fig. 14. — M 33, logio mass of third most massive cluster vs. logio SFR surface density. Symbols as in 
Figure [12] 
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Fig. 15. — M 33, logio mass of third most massive cluster vs. logio neutral gas surface density. Symbols as 
in Figure [12] 
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Table 5 

M 33. Fits to loGiq (Mard/MQ) vs. LOGio and a'^ 



PL 



Sgas/Mg pC 



3 62-63 (3.9 ±0.4) (0.6 ± 0.4) 

4 49-50 (3.5 ±0.6) (1.0 ±0.6) 

5 42-43 (4.0 ±0.8) (0.4 ± 0.8) 



3 62-63 (1.2 ±0.1) 

4 49-50 (1.1 ±0.2) 

5 42-43 (1.2 ± 0.2) 



(4 



Ssfr/Mq pc ^Gyr 



.18 ±0.03) 
(4.2 ±0.1) 
(4.1 ±0.1) 



3 62-63 (0.9 ±0.1) (3.8 ±0.1) 

4 49-50 (0.8 ±0.1) (3.8 ±0.1) 

5 42-43 (0.9 ±0.2) (3.7 ±0.2) 



NOTE.- 

each bin. 
(eq. 1). 



-Col. 
Col. 



(1): 
(3): 



number of bins, 
best-fit slope (eq 



Col. (2): number of clusters in 
1). Col. (4): best-fit intercept 



We examined two samples of clusters in M 51: 
one of 167 clusters 25 to 400 Myr old, and with a 
mass of at least IxlO^M©, and another one of 56 
clusters younger than 10 Myr, with masses of 5000 
Mq and up. For the older sample, regardless of the 
number of bins, wc find no correlation between the 
mass of the 3rd most massive cluster in each bin, 
Msrd, and the local surface densities of total gas, 
molecular gas, or star formation rate. There is, 
however, a hint of a correlation with local neutral 
gas surface density, such that logio M^^d — (0.4 ± 
0.2) logio Shi + (5.2 ± 0.3). 

For the younger sample, we find no correlation 
between Msrd, and either Eh2 or Ssfr- On the 
other hand, there is a correlation with Shi that is 
tighter and with a slightly steeper slope than the 
one found for old sample (although still consistent 
with it): logio M3,d = (0.6±0.1) logio Shi + (3.5± 
0.1). A correlation may also exist for the most 
recent burst of Msid with the total mass surface 
density Egas: logio A^Srd = (0.5 ±0.2) logio Sgas + 
(3.4 ± 0.2). If true, the comparison between the 
correlations with the older and younger burst im- 
plies that the normalization depends on the cluster 
mass and/or age ranges. 



The results could not be more different from 
the ones obtained for the flocculent galaxy M 33. 
From a 1-D radial analysis, we have previously de- 
rived for M 33: Afsrd oc ^ Msrd oc S 
]\ r vi-Oio.i 

J"3rd OC Zjc 



■^SFR 



^ y4.7±0.4. ^ yl.3±0.1. 

With a treatment that mimicks 



the 2-D analysis of M 51, we get for M 33: M^^d oc 

V3.8±0.3. M „ y^l.2±0.1. j./r ^ y0.9±0.1 

In M 33, the relation between M^j-d and SgpR 
is consistent with the expectations from pressure 
equilibrium considerations. On the other hand, 

the slope of the correlati on with Sgas is steeper 

than t he value assumed bv lPflamm-Altenburg fc Kroupa 



([20080 to explain the Ha cut-off in galaxy disks. 
There is no correlation with Shi, but M 33 is 
known to have been interacting, and the HI dis- 
tribution may not refiect the equilibrium configu- 
ration. The tight correlation with Eh2 and SsfRj 
and the looser one with Sgas agree with the recent 
results of a study of 33 nearby spiral galaxies by 



Schruba et al.l (|201ll ). who find a tight and lin- 
ear relation between Esfr and Ehs , with little 
dependence on Egas- 

How can one understand the diametrically op- 
posite results found for the grand-design spiral 
M 51? It is, of course, possible that stochastic 
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sampling is at work in M 51, such that the most 
massive cluster masses scale with the size of the 
sample. But then we would have to accept that 
stochastic sampling operates in some cases, like in 
M 51, and not in others, like in M 33. 

We propose a different hypothesis. The large 
azimuthal variations in gas and SFR surface den- 
sities in M 51 preclude a radial approach. Instead, 
we compared star cluster masses with local sur- 
face densities. The problem then is that, typi- 
cally, the clusters detected in the ACS data with 
masses above the completeness limit of IO^Mq are 
older than 10* yr (only 24 out of 167 have ages be- 
tween 40 and 90 Myr). Even if the cluster mass 
function is approximately independent of age over 
a timescale of a few hundred million years, clus- 
ters also complete between half and two full orbits 
around the galaxy in this time, and hence the lo- 
cal gas density around any given cluster probably 
has little to do with the gas density at the site and 
time of cluster formation. 

The two facts, that the correlation between 
Mard and Ehi seems to be tighter for the younger 
sample, and that this latter group of clusters might 
also exhibit a correlation between cluster mass and 
total gas surface density Sgas, give some credibil- 
ity to this hypothesis. Unfortunately, the resolu- 
tion of the radio data implies that we are com- 
paring clusters with sizes of the order of 10 pc 
with gas surface densities averaged over an area of 
w 750 X 750 pc^ . Even in the case of the younger 
sample, the measured Es are likely quite diluted 
compared to the actual densities relevant for the 
formation of the clusters. 

By contrast, 257 out of the 258 clusters kept 
for the final analysis in M 33 are younger than 
2.5x10^ yr (the age of the exception is 60 Myr) 
and, at any rate, the lack of a strong spiral pertur- 
bation means that the azimuthally averaged gas 
emission is in general representative of the star- 
forming conditions of clusters located at the same 
radii. 

To be able to perform a statistically signifi- 
cant analysis of clusters younger than 25 Myr, still 
close to their birthsites and whose local gas densi- 
ties resemble their star-forming environments, the 
depth and resolution of the cluster data need to be 
enough for the assembly of a sufficiently numer- 
ous sample above the mass function completeness 



limit F°l The gas data should also have a resolution 
comparable to the physical sizes of the clusters. 

The confirmation of the M 33 results will have 
to wait for the outcome of studies in closer targets, 
like the Magellanic Clouds, where the mass density 
distribution of clusters younger than a few x 10^ yr 
can be probed with large numbers down to 100- 
1000 Mq, and where the clusters' physical sizes 
are well matched by the resolution achievable at 
radio wavelengths; the data to do so already exist. 
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